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We construct a model for interacting electrons with strong spin orbit coupling in the pyrochlore 
iridates. We establish the importance of the direct hopping process between the Ir atoms and use 
the relative strength of the direct and indirect hopping as a generic tuning parameter to study the 
correlation effects across the iridates family. We predict novel quantum phase transitions between 
conventional and/or topologically non-trivial phases. At weak coupling, we find topological insulator 
and metallic phases. As one increases the interaction strength, various magnetic orders emerge. The 
novel topological Weyl semi-metal phase is found to be realized in these different orders, one of them 
being the all-in/all-out pattern. Our findings establish the possible magnetic ground states for the 
iridates and suggest the generic presence of the Weyl semi-metal phase in correlated magnetic 
insulators on the pyrochlore lattice. We discus the implications for existing and future experiments. 



I. INTRODUCTION 

Topological insulators^ ^ (TIs) have provided theorists 
and experimentalists alike with a new family of topo- 
logically non-trivial systems. In these materials, a suffi- 
ciently strong spin orbit coupling (SOC) leads to a pe- 
culiar band structure that cannot be adiabatically de- 
formed to that of a flat band insulator without closing 
the bulk gap. This leads to robust boundary states that 
display momentum-spin locking. The materials in which 
these gapless helical surface states have been observed are 
weakly interacting semiconductors, for which the above 
theory was constructed. An inviting question, therefore, 
relates to the kinds of quantum ground states that would 
arise in the presence of interactions in these systems or 
to interaction-driven TIs. For instance, several studies 
examined various kinds of fractionalized TIs.'^'^ 

In this context, transition metal oxides with 5d tran- 
sition metal elements may be ideal systems to search 
for TIs and new topological phases in the presence of 
interactions. In these systems, the strength of the in- 
teraction and that of the SOC are comparable, pro- 
viding a playground for the interplay between two ef- 
fects. In particular, the pyrochlore iridates, ^2^207, 
have been s uggeste d to host various topologically non- 
trivial states.S'^'^i'^ Here, A is a Lanthanide or Yit- 
trium, whose size affects the effective bandwidth of the 5c? 
electrons of Ir via the Ir-O-Ir bond angle, thereby tuning 
the effective strength of the interaction. Experiments on 
these compounds reveal metal-insulator transitions upon 
variation of temperature or chemicaP^ and exte rnal^ 
pressure, as well as indications of magnetism.!^^'^ 

In this work, we present a Hubbard-type model for the 
interacting electrons in the pyrochlore iridates and de- 
termine the ground state phase diagram using mean field 
and strong coupling methods. We find that it is impor- 
tant to include both the indirect hopping of 5d electrons 
of Ir through oxygens and the direct hopping between Ir 
sites. This is because the 5d orbitals of Ir are spatially 
extended and the nature of the ground state is sensitive 
to the relative strength of these hopping amplitudes. In 



the weakly interacting limit, both TIs and (semi-)metallic 
states are realized depending on the relative strength of 
different hopping amplitudes. This is in contrast to a 
previous work^ where only the indirect hopping process 
was considered and only the TI phase was obtained in the 
large SOC limit (with the ideal cubic crystal field). The 
interactions between electrons lead to two different mag- 
netically ordered ground states in different parameter re- 
gions. In particular, for inte rmediate interactions, the 
topological semi-metaPiElIIll (TSM) state with Weyl-like 
fermions appears in both kinds of AF phases. Our results 
suggest that the TSM state and the related Mott insulat- 
ing state can have different magnetic ordering patterns 
depending on the choice of the A-site ion or upon appli- 
cation of hydrostatic pressure, leading to the possibility 
of novel quantum phase transitions in the iridates. 



II. MODEL AND APPROACH 

In the atomic limit, the oxygen octahedra surrounding 
the Ir^+ ions create large cubic crystal fields that split 
the 5d orbitals into t2g and Cg multiplets. The five bd 
electrons of Ir**"*" occupy the t2g levels, leaving the high 
energy Cg levels empty. The angular momentum opera- 
tor projected into the t2g levels is effectively £ = 1 with 
an extra negative sign, i.e., Pt^^LPt^^ = —L^^^. The 
on-site SOC leads to a further splitting into an effective 
pseudospin j^s = 1/2 doublet and a jcs = 3/2 quadru- 
plet, the former lying higher in energy.!^ For sufficiently 
large SOC, the half-filled jcB = 1/2 doublets form a low 
energy manifold as the fully occupied jcS — 3/2 levels 
are sufficiently far from the Fermi level. In going to a 
tight-binding description, we need to take into account 
the different orientations of the local oxygen octahedra 
at each of the 4 sites in the unit cell (see Fig. |l(a)[ ). 
Previous studied^E^^^II considered nearest neighbour Ir-Ir 
hopping mediated by the oxygens. In this work, we also 
include the direct hopping between the Ir atoms, which 
is expected to be significant due to the large spatial ex- 
tent of Iridium's 5c? orbitals. We consider only the tt- 
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FIG. 1. a) Pyrochlore lattice of Ir atoms (large). The oxy- 
gens (small) bridging the Ir's in one unit cell are shown to- 
gether with the local axes they define, b) Phase diagram of 
the hopping Hamiltonian iifo- We set toxy = 1. The dashed 
line is = —21^/^. 



and cr-overlaps between the t2g orbitals, neglecting the 
usually smaller (5-overlap. This leaves us with two direct 
hopping parameters: t„ and It^. The resulting kinetic 
Hamiltonian reads 

Hq^ ^ {Tl]aa' +T^,aa')d'^Ria'^R'i'a' , (1) 

{RijR'i') ,aa' 

where R denotes the sites of the underlying Bravais FCC 
lattice of the pryrochlore lattice of Ir's, while i = I, . . . , 4 
labels the sites within the unit cell. The operator c?i?,i^(^) 
annihilates an electron in the pseudospin f (J,) state at 
site Ri. The two sets of matrices Tq and correspond to 
the oxygen mediatecP and direct hopping, respectively. 

We include interactions via an on-site Hubbard repul- 
sion between Iridium's d-electrons: 

H = Ho + Hu, (2) 
Hu ^ U^nRi^UR^i, (3) 

Ri 



where nma is the density of electrons occupying the 
bcff = 1/2, a) state at site Ri, with a =t, i- As we 
are interested in the magnetic phases expected at fi- 
nite [/, we perform a Hartree-Fock mean-field decou- 
pling Hu -J> -C^Ei?»(20'ffi) • jRi - Uri)^), where 

jRi = J2af!=t,i'^Ria'^'^0^R^0/'^ pseudospin oper- 

ator, whose expectation value will be determined self- 
consistently. We consider magnetic configurations pre- 
serving the unit cell so that (jRi) = (ji), i ~ 1,...,4, 
are the 4 order parameters under consideration. These 
are directly proportional to the local magnetic moment 
carried by the d-electrons. This follows from the fact 
that the projections of the spin and orbital angular mo- 
mentum operators onto the jeff = 1/2 manifold are pro- 
portional to the pseudospin operator: P^SP = — j/3 
and P^LP = -4j/3 with P = Pt,,Pi/2, where F*,^ 
projects onto the t2g subspace and P1/2 projects onto 
the Jeff = 1/2 subspace. This allows us to treat (ji) as 
the spontaneous local magnetic moment of the electrons. 



III. PHASE DIAGRAM 



A. Metal and topological insulator at U — 

We first examine the model at U = 0. Fig. 1 1(b)] shows 
the resulting phase diagram in terms of ta and (we set 
toxy = 1 throughout). Notice that both insulating and 
metallic phases exist. By virtue of the inversion symme- 
try of the crystal, we use the Fu-Kane formula^^ for the 
Z2 invariants in terms of the parity eigenvalues of the 
occupied states at the time reversal invariant momenta 
(TRIMs) to determine the topological class of each in- 
sulating phase. We find that both insulating phases are 
TIs with indices (1; 000). The TI phase adiabatically con- 
nected to io- = = corresponds to the large spin orbit 
limit of ReflHand is robust to the inclusion of weak direct 
hopping. As one tunes the direct hoppings, a metallic 
phase eventually appears by means of a gap closing at 
the r point. In the metal, the degeneracies at T become 
2-4-2 compared to 4-2-2 in the TI (with time-reversal 
and inversion symmetries all band are doubly degener- 
ate) . A similar situation occurs in Refs.T^ and f^TJ where 
a trigonal distortion of the oxygen octahedra drives the 
transition, not direct hopping as is the case here. The 
metallic phase is strictly speaking a semi-metal charac- 
terized by a point Fermi surface. Finite pockets can be 
generated by including very weak NNN hopping, as we 
have explicitly verified. Although we don't consider trig- 
onal distortions here, the direct hoppings alone can lead 
to qualitatively similar effects, e.g. the metallic phase 
resulting from the change in degeneracies at the F point. 
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B. Magnetic and topological phases at C/ > 



We now turn to the U > case. For convenience, 
we restrict our attention to a one-dimensional cut in 
the (t^jt^) space defined by = —2ta-/3, as shown in 
This is physically motivated since we expect 



Fig. EM 



and to have opposite signs, with the cr-overlap being 
the strongest. Moreover, the cut is representative as it 
intersects all the phases. In obtaining the finite U dia- 
gram, we performed an unconstrained analysis sampling 
over the space of all possible magnetic configurations pre- 
serving the unit cell. 

The resulting ground-state phase diagram appears in 
Fig. [2] First, we note that the TI is more resilient to the 
magnetic instability than the metal, as expected due to 
the presence of the bulk gap in the former. Second, the 
magnetic phase transition resulting from increasing U in 
the metal (TI) is second (first) order. Also, the magnetic 
order emerging from the TIs differs from the one found 
upon increasing U in the metal. In the latter case, we 
find an all-in/out configuration while in the former the 
ground state is 3-fold degenerate (modulo the trivial de- 
generacy j — >■ —j): all 3 states result from the all-in/out 
state by performing 7r/2-rotations on the moments in the 
unit cell. These rotations occur within either one of the 
planes bisecting the 3 triangles meeting at each corner of 
the tetrahedron. The order emergent in both TI states 
is the same. In section |IV[ we discuss how the differ- 
ent magnetic orders and the position of the transitions 
are actually connected to the corresponding ordering in 
the spin model obtained at large U: a.s ta is tuned, the 
induced Dzyaloshinski-Morya interaction alternates be- 
tween the only two symmetry allowed possibilities on the 
pyrochlore lattice, leading to different ordering. 



FIG. 3. Evolution of the spectrum as a function of U. At in- 
termediate U, in a), we can see a Weyl point along the T — L 
line, while in b), the spectrum naively seems insulating be- 
cause the Weyl points lie away from high symmetry fc-points. 
The dashed line is the Fermi level. 



C. Topological Semi- metal 

By examining the spectra of the ordered phases, we 
discover that the so-called topological semi-metal (TSM) 
is realized^'^ in the range io- > —1.67 and for a finite 
window of U. This semi-metallic phase has a Fermi "sur- 
face" composed of points, each with a linearly disper- 
sive spectrum of Weyl or two-component fermions, and 
may be considered as a 3D version of the Dirac points 
of graphene. The Hamiltonian near one such Weyl point 
takes the form 



H = Vq ■ q + ^Vi ■ qai 



(4) 



where q ~ k — ko is the deviation from the Weyl point 
at fep. The Pauli matrices ai represent the two bands 
involved in the touching, not (pseudo)spin. One can as- 
sign a chiral "charge" to these fermions, via the triple 
product of the 3 velocities: c — sgn(t)i • V2 x v^). The 
massless nature of the two-component Weyl fermions is 
robust against local perturbations, which is not the case 
in 2D. As explained in ReffTTi the only way to introduce 
a gap is to make two Weyl fermions with opposite chiral- 
ity meet at some point in the BZ. For this reason they 
are topological objects (see also the discussion below re- 
garding the surface states). Further details relating to 
the TSM can be found in Refs. [H [H [HI [H and [551 

The TSM appears in for both AF orders. In both cases 
we find a total of 8 Weyl points coming necessarily in 



4 




FIG. 4. (Colors online) Spectrum for a slab that is finite 
along the (110) direction, with ta = 0.25 and U — 2.8 which 
corresponds to a TSM. The inset shows the 2D BZ includ- 
ing the Fermi arcs connecting the projections of the Weyl 
points, where half of the Fermi arcs are located on the top 
surface, while the rest on the bottom one. The blue/red 
(dark/light) points correspond to Weyl bulk points with chiral 
charge +1/ — 1. 



a pair are made to move towards each other by increas- 
ing U, the corresponding Fermi arc shrinks, collapses to 
a point and disappears. 

In the TSM found at large io- , which we use to illustrate 
the Fermi arcs, there are no surface states along surfaces 
perpendicular to the (100), (010) or (001) directions. For 
these surfaces, the projection process onto the 2D BZ 
maps 3D Weyl points of opposite chirality onto the same 
2D fc-point. This leads to the absence of gapless surface 
states emanating from the 2D fc-point in question. For 
a surface perpendicular to the (110) direction, however, 
the projection is injective and Fermi arcs exist, as we 
illustrate in Fig. |4] 



4 inversion-symmetry related pairs. The location and 
migration of these Weyl points depends on the magnetic 
order. Let us first examine the TSM phase present in the 
all-in/out state. In this case, the 8 Weyl points are born 
out of the quadratic touching at the F point as the local 
moments spontaneously and continuously acquire a finite 
value with increasing U > Uc- Each pair of Weyl points 
lies on one of the four high symmetry lines joining T to 
the four L points, as can be seen in Fig. [3] For this reason 
we only get 8 touchings, in contrast to Ref.ll. where 24 
Weyl points are obtained. In their case they live off the 
high symmetry lines so that each point is tripled by the 
3-fold rotational symmetries about the F — L lines. Weyl 
points of opposite chirality annihilate at the 4 L points as 
U is increased. As they annihilate and create a gap, the 
parities of the highest occupied states at these TRIMs 
change sign. 

Let us now consider the TSM arising from the TI, 
where we again have 8 Weyl points. The major difference 
is that they do not occur along high symmetry lines, as 
can be seen in Fig. [3j We do not get 24 Weyl points 
because the magnetic order breaks the 3-fold rotational 
symmetries, which are preserved by the all-in/out state. 
We have explicitly located the Weyl points by looking at 
both the spectrum and density of states, which shows a 
characteristic {E — EpY scaling. The Weyl points don't 
annihilate at TRIMs, in contrast to the non-coUinear 
TSM. As a result there is no parity flip associated with 
the termination of the TSM phase when, upon increasing 
U , the system becomes insulating. 

Surface states: The non-trivial band topology of the 
TSM (each Weyl point is a monopole of the U(l) Berry 
connection) leads to chiral surface states on certain sur- 
faces, in analogy with the TI. In contrast with the lat- 
ter, the surface states of the former do not form closed 
Fermi surfaces, but rather open Fermi arcs. As argued in 
Refim the Fermi arcs join the projections of bulk Weyl 
points of opposite chirality. As bulk Weyl points forming 



IV. STRONG COUPLING EXPANSION 



In this section, we discuss the large U limit of our 
Hubbard Hamiltonian Eq. ([2]). We show how the effec- 
tive spin- 1/2 model obtained in that limit sheds light on 
the orders found in the mean field calculation as well as 
on the location of the phase transitions. In taking the 
limit where U is much larger than all hopping ampli- 
tudes (ioxy: io-j ^tt); wc cau usc second order perturbation 
theory to obtain the low energy spin Hamiltonian: 

H'^Y. ■ ^1 + ■iS^xS,) + 5f r°f 5^] (5) 

where the terms are, in order: the AF Heisenberg cou- 
pling, the Dzyaloshinski-Morya (DM) interaction and the 
anisotropic exchange. These correspond to the trace, an- 
tisymmetric and symmetric-traceless parts of the spin- 
spin interaction matrix, respectively. Let us focus on the 
bond between sites 1 and 2 (see Fig. |5(a)[ ), as the spin 
interactions for all other bonds can be determined using 
the crystal symmetries. We express the hopping Hamil- 
tonian between these two sites as 

Ht = -c\^ha,3C2f3 - cl^hlf^cip (6) 

where is a 2 by 2 complex matrix. Time-reversal 
symmetry restricts the matrix elements as follows: 

h = ta" + iv-(T (7) 

where t and v are real, and cr*^ is the identity matrix. 
We note that in order to derive the spin Hamiltonian, 
Eq. ([5]), we want to use the same quantization axes for 
both sites, i.e. we want the spin operators to be defined 
in the same coordinate system. 

Given the hopping matrix in the form Eq. ([t]), it can 
be shown quite simply that the Heisenberg, DM and 
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FIG. 5. &) D vectors corresponding to the "indirect" DM in- 
teraction. The "direct" type is obtained by changing the sign 
of all vectors, b) tVy as a function of the hopping strength t^. 
The sign of tVy determines the nature of the DM interaction: 
tVy > {tVy < 0) corresponds to the "indirect" ("direct") 
type. The red dots correspond to the metal-insulator transi- 
tions at ?7 = 0, which coincide with the points at which the 
nature of the magnetic order changes for U above the ordering 
threshold. 



anisotropic terms read 

rU 2 



D- = 2tv 

4 



U 



(8) 
(9) 
(10) 



If we turn to our microscopic hopping Hamiltonian, for 
the (1,2) bond we get 



t = a + btc 

V = Vy{0, 1, -1) 



(11) 



with Vy = a' + h't„ (12) 



where we have set toxy = 1 and = — 2i(j/3, as above. 
The coefficients a, 6, a', h' are positive rational numbers: 

a = 130/243 w 0.53 6 = 785/2916 ^0.27 (13) 
a' = 28/243 « 0.12 b' = 125/729 « 0.17 (14) 



We note that the v vector, hence D, is parallel to the 
opposite bond, (3,4), see Fig. |5(a)| This is a generic 
property of the pyrochlore lattice: as a consequence of 
crystal symmetry, a D vector for any given bond must 
be parallel to its opposite bond (in the sense that the 
4 sites form a tetrahedron) .I^Sl Moreover, if we know the 
DM vector for a single bond, crystal symmetries deter- 
mine the DM vectors for all other bonds in the lattice. 
Hence, there are only two possible sets of DM vectors 
{Dij}, called "direct" and "indirect". They are deter- 
mined by the sign of tVy for bond (1,2). (We could have 
picked another bond as the representative of the whole 
set.) The indirect (direct) type is defined as having the 
D vector for the bo nd bet ween sites 1 and 2 point along 
±(0, 1, —1). See Fig. 5 (a) [ for the configuration of D vec- 
tors corresponding to the indirect DM interaction. 

The nearest neighbour Heisenberg model together with 
a DM term on the pyrochlore lattice was studied by clas- 
sical Monte Carlo and mean field methods. •^^ First, the 
Monte Carlo study predicted a q = ordering, justify- 
ing the Ansatz used in the main text. Second, it was 
found that different magnetic orders arise depending on 
whether the DM interaction is of direct or indirect type. 
For the direct type, the configuration was found to be 
unique (up to time- reversal): the all- in/out order men- 
tionned above. Whereas for the indirect type, a continous 
manifold of degenerate orders was found, containing both 
coplanar and non-coplanar configurations. 

For the bond (1,2), we can extract from our micro- 
scopic Hamiltonian the value of the D vector: 



DA/U ^2tVy{S),l,-l) 



(15) 



Hence, \itvy = {a+bta-){a' +b'tcr) > we have an indirect 
exchange, otherwise it is direct. It is easy to see that the 
D vector changes direction when t = and Vy = 0, which 
correspond to t„ sa —1.99 and t„ ~ —0.67, respectively. 
For ta between these values, the DM interaction is of di- 
rect type, otherwise it is indirect. The behaviour of the 
DM interaction as a function of the direct hopping is 
shown in Fig |5(b)| We note that first value (-0.67) is al- 
most equal to the value at which the U = ground state 
goes from an insulator to a (semi)metal, tg- = —0.65. The 
magnetic orders we find for > —0.65 belong to the con- 
tinous manifold corresponding to the indirect DM term, 
while it is all-in/out when to- < —0.65, not too negative. 
Hence, the magnetic orders we get from our mean field 
calculation match those obtained in the strong coupling 
limit. The types of magnetic orders at intermediate U 
are found to be related to the type of DM interaction 
obtained in the large U spin model. 

We further note that the DM interaction becomes of 
indirect type for t„ < —1.99, which is sufficiently close to 
the second transition in the U — ground state, from the 
(semi)metal to the TI, which happens at ta = —1.67. For 
to- < —1.67, we get again the magnetic orders expected 
for an indirect DM interaction, again consistent with the 
large U limit. The bigger discrepency between the point 



6 



at which the D vector changes sign and the value of 
at which we observe a different ordering is probably due 
to the fact that the anisotropic exchanges increases in 
importance as t„ is increased, while it is smaller than the 
DM interaction near the first transition in the vicinity of 
t(j = 0. Hence, in that regime, we do not expect as good 
of an agreement with a spin model neglecting anisotropic 
exchange. 

V. DISCUSSION 

We have constructed a minimal (but sufficiently re- 
alistic) model to describe novel quantum ground states 
that may arise in the pyrochlore iridates. While not ap- 
preciated in previous works, it is shown that the inclu- 
sion of both indirect and direct hopping process of 5c? 
electrons of Ir is important in describing different mag- 
netically ordered states in the presence of interactions 
and their parent non-interacting ground states. A por- 
tion of our phase diagram is broadly consistent with a 
recent ah initio calculation,--' where upon increasing U , 
one encounters a metal, a topological semi-metal in the 
all- in/out magnetic configuration and finally a magnetic 
insulator. Since different choices of A-site ions in A2lr2 07 
lead to changes in both hopping amplitudes, our results 
suggest that different magnetic and topological ground 
states such as a topological insulator, the all- in/out and 



related AF states and various kinds of topological semi- 
metals, may arise in a variety of pyrochlore iridates. High 
pressure experiments on these compounds may reveal the 
intimate connection between the magnetic order in the 
stronger correlation regime and Tl/metal in the weak 
correlation limit, as theoretically explored in this work. 
For instance, recent transport measurements under high 
pressure!^ on Eu2lr207 indicate a continuous transition 
from an insulating ground state to a metallic one, mim- 
icking chemical pressure.'^ This could be connected to 
our continuous TSM-metal transition. Also, as the ex- 
istence of the TSM depends crucially on the magnetic 
order, it would be desirable to examine the effect of the 
magnetic fluctuations near the (semi-)metal-TSM tran- 
sition on thermodynamic and transport properties. 
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